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The mass flux of a binary mixture of smootii inelastic liard splieres or disks is determined by 
solving the Boltzmann equation by means of the Chapman-Enskog method to first order in the 
spatial gradients. As in the elastic case, the associated transport coefficients D, Dp and D' are given 
in terms of the solutions of a set of coupled linear integral equations which are approximately solved 
by considering the first and second Sonine approximations. The diffusion coefficients are explicitly 
obtained as functions of the coefficients of restitution and the parameters of the mixture (masses, 
diameters and concentration) and their expressions hold for an arbitrary number of dimensions. 
In order to check the accuracy of the second Sonine correction for highly inelastic collisions, the 
Boltzmann equation is also numerically solved by means of the direct simulation Monte Carlo 
(DSMC) method to determine the mutual diffusion coefficient D in some special situations (self- 
diffusion problem and tracer limit). As expected, the second Sonine approximation to D improves the 
predictions made from the first Sonine approximation. Beyond the above two limits, the quantitative 
variation of the complete set of diffusion coefficients across the parameter space is demonstrated 
by comparing the first and second Sonine approximations. A discussion on the convergence of the 
Sonine polynomial expansion is also carried out. Finally, as an application of the results derived 
here, the theoretical expressions for the diffusion coefficients are used to analyze segregation induced 
by a thermal gradient. The results obtained in this paper extend previous works carried out in the 
tracer limit (vanishing mole fraction of one of the species) by some of the authors of the present 
paper. 
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I. INTRODUCTION 



It is well established that granular matter under rapid 
flow conditions admits a hydrodynamic-like-description. 
At sufhciently low density, the Boltzmann kinetic equa- 
tion conveniently adapted to account for the inelastic 
character of collisions'"'^ has been used as the starting 
point to derive the corresponding hydrodynamic equa- 
tions. The essential assumption to get those equations is 
the existence of a normal solution^ defined to be one for 
which all the space and time dependence occurs through 
a functional dependence on the hydrodynamic fields. In 
the case of small spatial gradients, the Chapman-Enskog 
method^ provides a constructive means to get this normal 
solution and in particular, to obtain the Navier-Stokes 
(NS) constitutive equations in the first order of the ex- 
pansion. In this context, the study of hydrodynamics 
for granular gases follows similar steps as those made for 
ordinary gases. 

On the other hand, as in the elastic case,^ the ex- 
plicit form of the corresponding NS transport coefficients 
requires the solution of a set of linear integral equa- 
tions. The standard procedure of solving these integral 
equations consists of expanding the solutions in Sonine 



polynomials."* For simplicity, usually only the lowest So- 
nine polynomial (first Sonine approximation) is retained. 
However, in spite of this simple approximation, the re- 
sults obtained from this approach compare in general well 
with Monte Carlo simulations^"^ for mild degrees of in- 
elasticity. Although most of the theoretical results*"'^ 
have been devoted to monocomponent gases, some pro- 
gresses have been made in the past few years in the case of 
granular mixtures (namely, systems composed by grains 
of different masses, diameters and concentrations). In 
particular, in the context of granular mixtures at low 
density, Garzo and Dufty*'' have developed a kinetic the- 
ory which covers some aspects not completely covered 
in previous works. '"'"'^ Specifically, (i) the Garzo-Dufty 
theory goes beyond the weak dissipation limit so that 
it is expected to be applicable to a wide range of coef- 
ficients of restitution, and (ii) it takes into account the 
effects of nonequipartition of granular energy on the NS 
transport coefficients. As in the case of simple granular 
gases, the accuracy of the predictions of the Garzo-Dufty 
theory (which are based on the first Sonine approxima- 
tion) has been confirmed by numerical solutions of the 
(inelastic) Boltzmann equation by means of the direct 
simulation Monte Carlo (DSMC) method^° in the cases 
of the tracer diffusion coefficient^^'^^ and the shear vis- 
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cosity coefficient of a driven mixture. ^^'^"^ Discrepancies 
between theory and simulation appear to be important 
at strong dissipation for disparate mass and/or disparate 
size binary mixtures. Recently, the Garzo-Dufty theory 
has been extended to moderately dense binary mixtures^^ 
and the theoretical predictions compare also quite well 
with computer simulations. ^^'^^ 

A possible way of reducing the discrepancies between 
theory and DSMC results could be to consider higher- 
order terms in the Sonine polynomial expansion. In fact, 
recent works^^ analyzing diffusion of impurities in a gran- 
ular gas have shown that the second Sonine approxima- 
tion to the tracer diffusion coefficient yields a dramatic 
improvement (up to 50%) over the first Sonine approx- 
imation when impurities are lighter than the surround- 
ing gas in the range of large inelasticity. The results 
also show that the differences between the second Sonine 
approach and computer simulations are in general small 
(less than 4%) for arbitrarily large inelasticity. This good 
agreement stimulates the evaluation of the complete set 
of NS transport coefficients of a granular binary mixture 
(with arbitrary relative concentration) by retaining terms 
up to the second Sonine approximation. On the other 
hand, needless to say, the above goal is quite intricate 
due to the large number of collision integrals involved 
in the calculation. In this paper, we will cover partially 
this ambitious project by addressing the evaluation of the 
transport coefficients associated with the mass flux. 

We consider a binary mixture composed by smooth in- 
elastic disks {d = 2) or spheres {d = 3) of masses mi 
and 7712, and diameters cti and (T2- The inelasticity of 
collisions among all pairs is characterized by three inde- 
pendent constant coefficients of restitution an, a22, and 
ai2 = a2i, where < 1 is the coefficient of restitution 
for collisions between particles of species i and j. To first 
order in the spatial gradients, the constitutive equation 
for the mass flux (with i= 1,2) is given by^^ 

P p T 

(1) 

where D is the (mutual) diffusion coefficient. Dp is the 
pressure diffusion coefficient, and D' is the thermal dif- 
fusion coefficient. Here, n = ni + ^2 is the total num- 
ber density (rij is the number density of species i), 
p = mini -|- m2n2 is the total mass density, = rij/n 
is the concentration (or mole fraction) of species i, T is 
the granular temperature of the mixture and p = nT is 
the hydrostatic pressure. The main goal of this paper is 
to determine the diffusion coefficients D, Dp, and D' in 
terms of coefficients of restitution an, a22, and ai2 and 
the parameters of the mixture (relative masses, diame- 
ters and concentration). As said before and in contrast 
to our previous works, ^^'^^'^^ the diffusion coefficients will 
be explicitly obtained by considering contributions up to 
the second Sonine approximation. 

There are several reasons to address the above calcu- 
lation. First, given that the results reported in Refs. 
26 are limited to the tracer limit (xi — >■ 0), the ques- 



tion arises then as to whether (and if so, to what extent) 
the conclusions drawn before^ may be altered when 
one considers arbitrary concentrations. This goal is not 
only academic since, from a practical standpoint, many 
computer simulations^^'^^ and experiments^*^ in flowing 
granular mixtures involve finite conc;entrations. As a sec- 
ond reason, it must be noted that previous results ob- 
tained for ordinary mixtures^^ have clearly shown that 
while the first Sonine approximation can accurately de- 
scribe the shear viscosity and the thermal conductivity 
coefficients, it cannot achieve the same degree of accu- 
racy for the mutual and thermal diffusion coefficients. In 
this latter case, Kincaid et al.^^ concluded that the sec- 
ond Sonine approximation is much better approximation 
than the first one for a wide range of values of masses 
and sizes. A third motivation to improve the evalua- 
tion of the NS transport coefficients lies in the fact that 
the reference homogeneous cooling state (HCS) is known 
to be unstable against long wavelength spatial pertur- 
bations, leading to vortex and cluster formation. Since 
this instability can be well characterized"^^ through a lin- 
ear stability analysis of the hydrodynamic equations, a 
more accurate evaluation of the NS transport coefficients 
for large inelasticity may help to understand the physi- 
cal mechanisms involved in this instability. Finally, as a 
fourth motivation and given that the second Sonine ap- 
proach is expected to differ from the first one at strong 
dissipation, the results reported here can be of practical 
interest since the range of high inelasticities has growing 
interest in experimental works^* and is also exhibited by 
wetted particles,''''' which are critical in processes such as 
coagulation, filtration or pneumatic conveying. 

The plan of the paper is as follows. First, in Sec. II 
The Boltzmann equation and its corresponding balance 
hydrodynamic equations for the mass, momentum and 
energy are recalled. In Sec. Ill, the diffusion transport co- 
efficients D, Dp, and D' are given in terms of the solution 
of a set of coupled linear integral equations previously de- 
rived by Garzo and Dufty.^^ These integral equations are 
approximately solved by using the first and second So- 
nine approximations, where explicit forms for the above 
transport coefficients are provided. Technical details of 
the calculations carried out in this paper are relegated to 
two Appendices. Next, the theoretical approaches (first 
and second Sonine approximations) are compared in Sec. 
IV with available and new simulation data obtained from 
numerical solutions of the Boltzmann equation by using 
the DSMC method for the self-diffusion and tracer dif- 
fusion coefficients. Two- and three-dimensional systems 
are considered. The dependence of the complete set of 
diffusion coefficients on the parameter space of the sys- 
tem is widely analyzed in Sec. V. The results indicate 
that, while the convergence of the Sonine expansion on 
the mutual diffusion coefficient D is relatively good, in 
the cases of the coefficients Dp and D' is not as good. 
The expressions derived here for the diffusion transport 
coefficients are applied to analyze thermal diffusion seg- 
regation in Sec. VI. This is a quite interesting problem 
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not only from a fundamental point of view but also from a 
more practical point of view. The influence of the second 
Sonine solution on the phase diagrams for segregation is 
assessed in several cases showing that previous results^^ 
obtained by using the simple first Sonine correction are 
significantly affected. We close the paper in Sec. VII with 
a brief discussion of the results. 



species, which mcasiire their mean kinetic energies. The 
exact macroscopic balance equations for ni{r,t), u(r,t), 
and T{r,t) follow directly from Eq. (2) (and its corre- 
sponding counterpart for /2) by multiplying with 1, mjV, 
and ^niiV^ and integrating over v. They are given by^^ 

DtUi + riiV • u + = , (7) 

rui 



II. BOLTZMANN KINETIC THEORY FOR 

GRANULAR BINARY MIXTURES Au + p" V • P = , (8) 



Let us consider a granular binary mixture where the 
density of each species is sufficiently low. In this case, 
all the rck^vant information on the state of the mixture 
is accounted for by the velocity distribution functions 
/i(r, v;t) of each species {i = 1,2). These distributions 
are accurately described by the set of nonlinear Boltz- 
mann equations^ 

(at + v-V)/i(r,v,t) = Jii[v|/i(i),/i(<)] 

+ Jl2[v|/l(i),/2(t)] , (2) 

and a analogous equation for /2(r, v; t). The Boltzmann 
collision operators [v|/i, fj] are given by 

Jij [Vll/i, fj] = erf-' J dV2 J d^ • gl2)(^ • gl2) 



X K//^(r,v;,t)/j(r,v^,t) 
-/j(r, vi,i)/,(r,V2,t)] , 



(3) 



where d is the dimensionality of the system, = 
((Ti + (jj) /2, is a unit vector along the line of centers, 
Q is the Heaviside step function, and gi2 = Vi — V2 is 
the relative velocity. The primes on the velocities denote 
the initial values {v'^, Vj} that lead to {vi,V2} following 
a binary (restituting) collision: 



V'l = Vi - Hji (1 + a^/) (ct • gi2)CT, 

V2 = V2 + iiij (1 + a:r.') {a ■ gi2)CT, 



(4) 



where fiij =mi/ {rrn +mj). 

In the case of granular mixtures, the relevant hydro- 
dynamic fields are the number densities ni{r,t), the flow 
velocity u(r, t), and the granular temperature T(r, t). In 
terms of the velocity distribution functions /i(r, v, t), the 
above fields are defined respectively as 



ni= dv/i(v) , pu 



2 

E 



nii 



dvv/i(v) , (5) 



^ = E-.^. = E^/dvl^V.(v), (6) 

i=l i=l 

where V = v— u is the peculiar velocity. The third equal- 
ity of Eq. (6) defines the kinetic temperatures Tj for each 



A^-?E^ + i(V-q+P:Vu) = -CT. (9) 



n rrii 



In the above equations, Dt = 9t + u • V is the material 
derivative, 

ji=mij dvV/i(v) (10) 
is the mass fiux for species i relative to the local fiow. 



2 

E 

i=l 



rrii 



dvVV/i(v) 



(11) 



is the total pressure tensor, 

2 



Em f 



dvy2v/,(v) 



(12) 



is the total heat fiux, and 



C = Y: xmCi = EE ^ / dvF^^[v|/., /.] , 
i=l p i=l j=l " 

(13) 

is the total "cooling rate" due to inelastic collisions 
among all species. In Eq. (13), 7i = Ti/T and the sec- 
ond equality defines the "cooling rates" Q for the partial 
temperatures Tj. 

The balance equations (7)-(9) do not constitute a 
closed set of equations for the hydrodynamic fields unless 
one knows the functional dependence of jj, P, q, and C 
on the above fields. On the other hand, for times longer 
than the mean free time, the distribution functions /j 
are expected to adopt the form of a normal or hydrody- 
namic solution such that all space and time dependence 
of fi occurs through the hydrodynamic fields: 

fi{r,v,t) = f[y\xi{r,t),p{r,t),T{r,t),u{r,t)] . (14) 

Note that we have taken the set {xi,p, T, u} as the d + 3 
independent fields of the two-component mixture. As 
mentioned in Ref. 27, in the case of inelastic systems, 
there is more flexibility than in ordinary mixtures to 
chose the set of relevant hydrodynamic fields since the 
specific set of gradients contributing to each fiux is only 
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restricted by fluid symmetry considerations. Here, as in 
our previous works for dilute granular mixtures, -'^^'^^ we 
have chosen the set {a;i,p, T, u} since they are the most 
accessible fields from an experimental point of view. In 
particular, a contribution proportional to Vp (which is 
absent in the elastic case) appears in the mass and heat 
fluxes. 

In the case of small spatial variations (i.e., low Knudsen 
numbers), the functional dependence (14) can be made 
local in space through an expansion in the gradients of 
the hydrodynamic fields. This is the procedure followed 
in the Chapman-Enskog method^ to get an approximate 
solution to the Boltzmann equation. Thus, the distribu- 
tions /, are written as 



(1) 



f(2) 



(15) 



where each factor of e (formal nonuniformity parame- 
ter) means an implicit gradient of a hydrodynamic field. 
In the first-order of the expansion, the NS constitutive 
equations for the mass, momentum and heat fluxes can 
be derived. In this paper, we will focus on our attention 
to the first-order contribution j -^^ to the mass fiux. 



III. DIFFUSION TRANSPORT COEFFICIENTS 

The application of the Chapman-Enskog method to the 
Boltzmann equation allows one; to dc!tc;rminc the form of 
the NS transport coefficients of the mixture. In particu- 
lar, the mass flux j-^' is given by Eq. (1) where the dif- 
fusion transport coefficients D, Dp, and D' are defined, 
respectively, as 



D = -- 



dm2n 



I 



dvV-Au 



Dp = -!^ / dvV -Bi, 



(16) 



(17) 



_^(0) ^TdT+pdp)+Ci-2C^°^ 



P 

_^(o) (^Tdr + pdp) + £2 - 2C^°^ 
H 1'2, 



2T 



-C(°) {TdT +pdp)+ £2 -IC^°^ 
2T 



02, 



Bi+A^iB2 = Bi 

(21) 

B2+M2B1 = B2 
(22) 

ei + Mie2 = Ci 

(23) 

C2+M2C1 = C2 

(24) 



Here, ('^^^ = ([^'' = (^If^ is the cooling rate evaluated with 
the zeroth-order distribution and we have introduced the 
linearized Boltzmann collision operators 

C,X = - (Jn[/f + Jn[X, + J,2[X, /(°^]) , 

(25) 



MiX = -J,2[fi°\x]. 



(26) 



The corresponding forms for the operators £2 and M2 
can be easily obtained from Eqs. (25) and (26), respec- 
tively, by just making the changes lo2. In addition. 



d 



B.(V) 



nT 



f(0) 



(27) 

(28) 



D' 



miT 
dp 



I 



dvV Ci. 



(18) 



As in the case of elastic collisions,^'^^ the quantities Ai, 
Bi, and d {i = 1,2) are the solutions of the following set 
of coupled linear integral equations: "'^^ 



_^(o) (TdT+pdp)+£i 



A\ -h A^i^2 = Ai 



-C(°^ {TdT + pdp) + £2] A2 + M2AX = A: 



(19) 



(20) 



Ci(V) 



1 

T 



2dN 



(v/(o)) 



V. 



(29) 



It is worthwhile remarking that so far the expressions 
for the transport coefficients D, Dp, and D' are exact. 
However, in order to determine the dependence of the 
above coefficients on the parameters of the mixture, one 
needs to solve the integral equations (19)-(24) and to 
know the explicit form of the (local) HCS distributions 
fj;^\ With respect to this latter point, both theoretical"^*^ 
and computer simulation^^'^^ results have shown that 

in the region of thermal velocities //"^(V) is well rep- 
resented by its Maxwellian form 



/f (V)-^Am(V) 



d/2 



exp 



2Ti 



(30) 
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Thus, in order to get simple and accurate expressions for 
the diffusion transport coefficients, we will neglect here 

the non-Gaussian corrections to /i'^^(V). While these 

corrections arc not important in the case of the mass flux 
and the pressure tensor, ■^^ the impact of them on the heat 
flux is not negligible in highly dissipative gases/'^"'^ Ac- 
cordingly, a theory incorporating the above non-Gaussian 
corrections does not seem in practice necessary for com- 
puting the diffusion transport coefficients. 

Regarding the unknowns Ai, Bi, and Cj, the standard 
method consists of approximating them by Maxwellians 
(at different temperatures) times truncated Sonine poly- 
nomial expansions. For simplicity, usually only the low- 
est Sonine polynomial (first Sonine approximation) is 
rctaincd^'^^'^"^ and the results obtained from this simple 
approach agrees in general relatively well with numeri- 
cal results^''''''' for granular mixtures obtained from the 
DSMC method. However, as for ordinary mixtures, "^^ 
significant discrepancies between theory and simulation 
appear when one considers disparate values of mass and 
diameter ratios at small values of the coefficients of resti- 
tution. This disagreement could be mitigated in part if 
one considers higher-order terms in the Sonine polyno- 
mial expansion, as has been reported in detail^^'^® in the 
tracer limit {xi 0) at the level of the diffusion coeffi- 
cient D. In particular, as said in the Introduction, it is 
shown that the accuracy of the second Sonine approxima- 
tion for D is much better than the first Sonine approxi- 
mation when the tracer particles are lighter than the par- 
ticles of the gas. Motivated by these results, our goal here 
is to civaluatc the complete set of diffusion coefficients D, 
Dp, and D' up to the second Sonine approximation as 
functions of the coefficients of restitution (an, 0:22, and 
Q12) and the parameters of the mixture (masses m^, di- 
ameters cTj and concentration xi). Therefore, the present 
analysis generalizes to arbitrary concentration our previ- 
ous theoretical results derived in the simple tracer limit 
case. 

In the second Sonine approximation, the quantities 
Bi, and Ci are approximated by 

Ai(V) ^ /i,M [ai,iV + ai,2Si(V)] , 

A2(V) ^ /2,M [a2,lV + a2,2S2(V)] , 

Bi(V) ^ /i,M [61,1V + 61,281 (V)] , 

B2(V) ^ /2,M [62,1V + 62,2S2(V)] , 
Cl(V) ^ /i,M [C1,1V + Ci,2Si(V)] , 



where 



C2(V) ^ /2,M [C2,lV + C2,2S2(V)] , 



s,(v)=('im.y2-^r.)v. 



(31) 

(32) 
(33) 
(34) 
(35) 
(36) 

(37) 



The coefficients {a^^i, Ci,i} are related to the trans- 
port coefficients D, Dp, and D' , respectively, as 



n2T2 mim2n ^ 
ai,i = 7^0-2,1 = T^D, 



61,1 = -—^62,1 = - 



pniTi 



(38) 
(39) 

niTr^'^ TniTi" " ^^'^^ 
Upon writing the first equalities in Eqs. (38) (40) use has 
been made of the property j^^^ = — j2^^ The coefficients 
{(ii,2,bi^2,Ci,2} are defined as 



ci,i = — 



n2?2 



-C2,i = 



P 



-D'. 




d{d + 2) mTl 



(ivSi(V) 



(41) 



The diffusion transport coefficients D, Dp, and D' and 
the second Sonine coefficients 0^,2, 6j,2, and Ci,2 are deter- 
mined by substitution of Eqs. (31)-(36) into the integral 
equations (19)-(24), multiplication by m^V and Si(V), 
and integration over velocity. The procedure is lengthy 
and follows similar mathematical steps as those made 
before^^'^^ in the tracer limit {xi — >■ 0). Technical details 
on this calculation have been relegated to the Appendix 
A. 

For the sake of convenience, we introduce dimension- 
less forms for the diffusion coefficients as 



D = 



pT 



-D* 



mim2fo 



r?T 

D =—D* 
pfo ^ 



D' 



nT 

PVQ 



D' 



where 



vo = \/7rn<Ti^\/2T 



^mi -I- m2 
mim2 



(42) 



(43) 



is an effective collision frequency. The expressions of the 
(reduced) second Sonine approximations £'*[2], -D*[2], 
and D * [2] are obtained from Eq. (A34) where one has to 
invert the 9x9 square matrix H defined by Eq. (A29). 
The solution gives the (reduced) diffusion coefficients 
for any dimension d as a function of the mole frac- 
tion x\, the mass ratio jj, = mi/ 1112, the diameter ra- 
tio uj = cri/<T2, and the coefficients of restitution an, 
a22, and ai2 = a2i. The explicit forms of D*[2], £>*[2], 

and D * [2] are too large to be displayed here and will be 
omitted for the sake of simplicity. On the other hand, 
since j^^^ = -j2"^\ D*[2] must be symmetric while D*[2] 
and D'*[2] must be antisymmetric with respect to the 
exchange 1 o 2. We have checked that our expressions 
verify the above symmetry properties. 

It is quite apparent that the influence of the parameters 
of the mixture on the second Sonine approximations is 
rather complicated, given the large number of parameters 
involved in the system. Thus, in order to show more 
clearly the dependence on each parameter on diffusion, 
it is instructive to consider first some simple cases. 
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A. Some special limits 

Let us first consider the first Sonine approximations 

D*[l], D*[l], and D'*[1]. They can be obtained from the 
general results by taking aj^2 = ^i,2 = Cj,2 = 0. In this 
case, one gets 



D*[l] = 



d 



+ 



dC 
dxi 



D;[1] = XI 71 



X2 + jJiXi 



2v* I ^ 



-c + — 



(44) 



(45) 



Dp through Eq. (38)) and cn (which defines the thermal 
diffusion coefficient D' through Eq. (38)) can be given 
in terms of 622 and C22. The corresponding expressions 
for Z?*^], -D*[2], and D'*[2] coincide with those previ- 
ously obtained'^^ by following an independent route. In 
particular, the tracer diffusion coefficient £)*[2] is given 

by 



D*[2]=7 



(z^*-K*)Ki-ic*)-'^iVri' 



(49) 



where u* is given by Eq. (Bl) with Xi = and the quan- 
tities z/j^j, and rj^i are defined in the Appendices A 
and B. 

All the above limits confirm the self-consistency of the 
results derived in this paper for the second Sonine ap- 
proximation to the diffusion coefficients D, Dp, and D' 
of a granular binary mixture. 



(46) 



where C* = C^^V'^o and u* is given by Eq. (Bl). The 
temperature ratio 7 = T1/T2 is determined from the con- 
dition CI = C2 = C*: where the partial cooling rates Q 
are given by Eq. (B13). The expressions (44)-(46) agree 
with those derived in previous works. ^"^'^^ 

Another interesting situation is the case of mechani- 
cally equivalent particles (mi = m2, (Ti = (T2. oii = 
a22 = 0:12 = a). In this simple situation, as expected, 
our results yield D*p[2] = D'*[2] = and 



D*[2] = D*[l] 



1 + a 
d 



12a'^ + 3{2d - 3)a + 8 + lOd 
^ 12a3 + {6d - 5)a2 + {16d + l)a + lOd + 12 ' 

(47) 

where the first Sonine approximation £'*[1] is simply 

2r(i) d 



D*[l] = 



(l + a)2 



(48) 



As expected, the expression of the self-diffusion coeffi- 
cient D* [2] holds for any relative number of tagged parti- 
cles since it is independent of xi. Equation (48) coincides 
with previous results for the self-diffusion coefficient. 

Let us consider finally the tracer limit, namely, wc as- 
sume that the concentration of one of the species (say for 
instance, species 1) is negligible {xi — >■ 0). In this limit, a 
careful analysis of the matrix equation (A28) defining the 
Sonine coefficients Uij, bij, and Cij shows that 022 = and 
the coefficients an (which defines the diffusion coefficient 
D through Eq. (38)) and ai2 arc decoupled from the re- 
maining 6 Sonine coefficients. Moreover, the coefficients 
622 and C22 associated with the excess component also 
verify an autonomous set of equations so that, the coeffi- 
cients 611 (which defines the pressure diffusion coefficient 



IV. COMPARISON WITH DSMC RESULTS 

Needless to say, the improvement of the analytical re- 
sults by considering the second Sonine approximation 
for the diffusion coefficients is not completely guaran- 
teed unless the Sonine polynomial expansion is conver- 
gent. The analysis of the convergence of the Sonine poly- 
nomial expansion for granular gases is a quite intricate 
problem due to the mathematical difficulties of obtaining 
higher-order Sonine corrections to the transport coeffi- 
cients. Thus, the works devoted to this issue in granular 
systems arc more scarce than for ordinary gases where 
for instance the analysis of transport properties of dense 
binary mixtures**^ with one tracer component {xi 0) 
indicates that the convergence of the Sonine expansion 
improves with increasing values of the mass ratio mi/m2. 
In this Section, we will compare the first and second So- 
nine approximations of the mutual diffusion coefficient 
D with computer simulation results obtained by numer- 
ically solving the Boltzmann equation by means of the 
DSMC method. As in previous studies,^^'^^'^^ due to 
the difficulties for measuring the coefficient D for gen- 
eral values of the mass ratio and the mole fraction, we 
will consider the self-diffusion (mi = m2) and tracer dif- 
fusion {xi — )• 0) coefficients. Moreover, to cover more 
general systems than those considered in our previous 
simulations, -^^'-^^'^^ we will assume that q;i2 ^ q;22 when 
the intruder and the gas particles are mechanically dif- 
ferent. 

The adaption of DSMC method to analyze binary 
granular mixtures has been described in previous works 
(see, for instance, Ref. 28), so that here we shall only 
mention some aspects related to the diffusion of impu- 
rities in a granular gas under HCS. In the tracer limit 
(ni <C 712), during our simulations collisions 1-1 are not 
considered, and when a collision 1-2 takes place, the post- 
coUisional velocity obtained from the scattering rule is 
only assigned to the tracer particle (species 1). Accord- 
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FIG. 1: Plot of the (reduced) self-diflfusion coefficient D{a)/D{l) as a function of the coefficient of restitution a as given by the 
first Sonine approximation (dashed Une), the second Sonine approximation (sohd Une), and Monte Carlo simulations (symbols). 

Here, -D(l) is the clastic value of the self-diffusion coefficient consistently obtained in each approximation. The left panel is for 
hard disks (d = 2) while the right panel is for hard spheres (d = 3). 



ing to this scheme, the numbers of particles have simply 
a statistical meaning and can be arbitrarily chosen. 

In an unbounded system like ours, the DSMC method 
has two steps that are repeated in each time iteration. 
In this first step, the system (tracer and gas particles) 
evolves from the initial state to the HCS. In the second 
step, the system is assumed to be in the HCS and then 
the diffusion coefficient D(t) is measured from the mean 
square displacement of the impurity as 

(50) 

Here, |r(f) — r(0)| is the distance traveled by the impurity 
from t = imtil time t, t = being the beginning of 
the second step. Moreover, (• • • ) denotes the average 
over the N impurities and St is the time step. In our 
simulations, we have typically taken a time step St = 
2.5 X 10~'^iy~^ and N = 2 x 10^ simulate d particles for 
each species. Here, v = n2a'^~^ \/2T/m2 is an effective 
collision frequency for gas particles. 

We will consider first the self-diffusion coefficient, 
which is independent of the mole fraction xi [sec Eqs. 
(47) and (48)]. The simulation data obtained from 
DSMC method along with both Sonine approximations 
for the reduced coefficient D{a)/ D{1) are presented in 
Fig. 1 for disks {d = 2) and spheres {d = 3). Here, D{1) 
refers to the elastic value of the self-diffusion coefficient 
consistently obtained in each Sonine approximation. The 
data corresponding to d = 3 for a > 0.5 and d = 2 for 
a > 0.6 were reported in Refs. 21 and 22, respectively, 
while those corresponding to d = 3 and d = 2 for a < 0.5 
have been obtained in this work. It is quite apparent 
that the first Sonine approximation performs well for not 
strong values of dissipation, but the agreement between 
theory and simulation improves over the complete range 
of values of the coefficient of restitution when the second 
Sonine approximation is considered (especially for hard 
disks). This confirms again the accuracy of the second 
Sonine approach even for quite extreme values of dissi- 
pation. 



1.2 F 




ai2 

FIG. 2: Plot of the (reduced) mutual diffusion coefficient 

D{ai2)/D{l) versus the coefficient of restitution ai2 in the 
tracer limit {xi — )• 0) for a granular gas of hard spheres with 
Ld — 1/2, /i = 1/4 and 022 = 0.5. The dashed and solid 
lines are first and second Sonine approximations, respectively, 
while the symbols are the Monte Carlo simulation results. 
Here, -D(l) is the elastic value of the mutual diffusion coeffi- 
cient consistently obtained in each approximation. 



Consider now the situation in which impurities and 

particles of the gas are mechanically different (i.e., they 
can differ in size, mass and coefficients of restitution). 
Although not shown here, as expected, comparison be- 
tween theory and simulation shows that the Sonine poly- 
nomial expansion exhibits a better convergence (namely, 
both Sonine approximations compare well with numeri- 
cal results) when the impurity is heavier and/or larger 
than the gas particles while this convergence is worsen as 
fi and/or 00 significantly decreases. These findings agree 
with the conclusions obtained for elastic collisions.^" To 
illustrate this behavior. Fig. 2 shows the dependence 
of the ratio D{ai2)/ D{1) on the coefficient of restitu- 
tion ai2 for hard spheres with w = 1/2, fj, = 1/4 and 
0:22 = 0.5. The present comparison complements previ- 
ous results^^'^^'^^ reported for the special case ai2 = a22- 
We observe that the first Sonine approximation clearly 
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FIG. 3: Plot of the (reduced) mutual difFusion coefficient 
D/D{1) as a function of the mass ratio in the tracer limit 
{xi — >■ 0) for a granular gas of hard spheres with tj = 1/2 and 
a (common) coefficient of restitution a = 022 ~ Q12 = 0.5. 
The dashed and solid lines are first and second Sonine ap- 
proximations, respectively, while the symbols are the Monte 
Carlo simulation results. Here, D{1) is the elastic value of 
the mutual diffusion coefficient consistently obtained in each 
approximation. 



us to assess the reliability of the second Sonine solution 
to D beyond the tracer limit. The fact that the sec- 
ond Sonine expression for D in the self-diffusion problem 
(which holds for any value of xi) compares quite well 
with DSMC results suggests that the good agreement 
found for xi — > would be also kept for arbitrary values 
of the mole fraction, even when both species are mechan- 
ically different. More simulations are needed to support 
the above expectation. 
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overestimates the simulation results while the second So- 
nine approximation to D(ai2) exhibits good agreement. 
On the other hand, the quantitative discrepancies be- 
tween the second Sonine solution and simulation data 
are larger than those observed for the self-diffusion prob- 
lem (see Fig. 1), especially for strong dissipation. Thus, 
one perhaps would have to consider the third Sonine cor- 
rection to obtain a better prediction for the diffusion co- 
efficient. 

We explore now the influence of the mass ratio /i on 
the accuracy of the two first Sonine approximations. Fig- 
ure 3 shows the ratio D/D{1) versus the mass ratio for 
hard spheres with w = 1/2 and a (common) coefficient 
of restitution a = a22 = c«i2 = 0.5. We find that the 
second Sonine approximation D[2] differs form the first 
Sonine approximation £'[1] as the mass ratio fj, is varied. 
For the system studied in Fig. 3, the disagreement be- 
tween both approaches turns out to be significant when 
the impurity is heavier than the gas particles. Thus, 
for instance when /i = 5, the first Sonine approxima- 
tion to the ratio D/D{1) differs by 26% from the sec- 
ond Sonine approximation. The comparison with simu- 
lation data shows again that the theoretical predictions 
are clearly improved when one takes the second Sonine 
solution. However, the quantitative differences between 
the second Sonine solution and DSMC results seem to in- 
crease as the mass ratio increases. In this case, as in Fig. 
2, one should consider higher-order terms in the Sonine 
polynomial expansion to get a more accurate approach. 

The results reported in this Section confirm again the 
reliability of the second Sonine approximation for the 
mutual diffusion coefficient D, at least in the cases of 
self-diffusion and tracer limit. Unfortunately, the lack of 
available simulation data for finite mole fraction prevent 



FIG. 4: (color online) Plots of the reduced coefficient 
D(q)/D(1) as a function of the (common) coefficient of resti- 
tution a for hard spheres with xi = 0.2, ai = (T2 and two 
different values of the mass ratio fi = mi/m2. The solid lines 
correspond to the results obtained from the second Sonine 
approximation, the dashed lines refer to the (standard) first 
Sonine approximation and the dotted lines correspond to the 
modified first Sonine approximation. Here, D{1) is the elastic 
value of D consistently obtained in each approximation. 
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FIG. 5: (color oidine) Plot of the reduced coefficient 
Dp{a)/ Dp{l) as a function of the (common) coefficient of 
restitution a for hard spheres with xi — 0.2, ai — 02 and 
two different values of the mass ratio ^ = m\/m2. The solid 
lines correspond to the results obtained from the second So- 
nine approximation, the dashed lines refer to the (standard) 
first Sonine approximation and the dotted lines correspond to 
the modified first Sonine approximation. Here, -Dp(l) is the 
elastic value of Dp consistently obtained in each approxima- 
tion. 
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FIG. 6: (color online) Plot of the reduced coefficient D * (a) as 
a function of the (common) coefficient of restitution a for hard 
spheres with x\ = 0.2, ai = (J2 and two different values of the 
mass ratio /i = m\/m2. The solid lines correspond to the 
results obtained from the second Sonine approximation, the 
dashed lines refer to the (standard) first Sonine approximation 
and the dotted lines correspond to the modified first Sonine 
approximation. 



FIG. 9: (color online)The ratio of the second and first So- 
nine approximations D'[2]/Z)'[l] to the thermal diffusion co- 
efficient versus the mole fraction x\ for uj = 1, a = 0.8 and 
two values of the mass ratio = A and ^ = 1/3). 



V. DEPENDENCE OF THE DIFFUSION 
COEFFICIENTS ON THE PARAMETERS OF 
THE MIXTURE 




FIG. 7: (color online)The ratio of the second and first Sonine 
approximations 7?[2]/D[l] to the mutual diffusion coefficient 
versus the mole fraction xi for u} = 1, a = 0.8 and two values 
of the mass ratio (/.i = 4 and /i = 1/3). 
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FIG. 8: (color online)The ratio of the second and first Sonine 
approximations Dp [2] /Dp [1] to the pressure diffusion coeffi- 
cient versus the mole fraction xi for u} = 1, a = 0.8 and two 
values of the mass ratio (/x = 4 and fj, — 1/3). 



Once the reliability of the second Sonine solution to 
the mutual diffusion coefficient D has been confirmed 
in the previous Section, our goal now is to provide a 
systematic study of the dependence of the complete set 
of diffusion coefficients D, Dp and D' on the parameter 
space of the system. However, as mentioned before, the 
first and second Sonine approximations to the (reduced) 
diffusion transport coefficients depend on many parame- 
ters: {xi, TOi/to2, cri/cr2, ail, q;22, 0112}. This complexity 
also occurs for ordinary binary mixtures, except that now 
we have the coefficients of restitution as additional pa- 
rameters. Also, to reduce the number of independent 
parameters, the simplest case of a common coefficient of 
restitution (an — a22 = ai2 = a) and a common di- 
ameter (cTi — (72 ) is considered. The latter assumption 
is justified because the dependence of D*, D* and D * 
on the diameter ratio uj is very weak. Moreover, hence- 
forth we only analyze the physical case of hard spheres 
(d = 3) and so, the parameter space is reduced to three 
quantities: {xi, mi/m2, a}. 

The first and second Sonine approximations of the (re- 
duced) transport coefficients D{a)/D{l), Dp{a)/Dp{l), 
and D *(a) are plotted in Figs. 4, 5, and 6, respectively, 
for xi = 0.2 and two values of the mass ratio /i. The 
diffusion coefficients have been reduced with respect to 
their elastic values (consistently obtained in each Sonine 
approximation), except the thermal diffusion coefficient 
D' since it vanishes for elastic collisions when one consid- 
ers the first Sonine approximation. In this latter case, we 
have plotted the reduced coefficient D * defined by the 
third relation in Eq. (42). For the sake of comparison, 
we have also included the results derived from a mod- 
ified version of the first Sonine approximation.^^ This 
approach consists of replacing the Maxwellian distribu- 
tion in the first Sonine solution by the HCS distribution. 
Figure 4 shows the a-dependence of the mutual diffusion 
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coefficient obtained from tlie tlrree different approxima- 
tions (standard and modified first Sonine approximation 
and tlie second Sonine approximation) for two mass ra- 
tios. We observe that the first Sonine approximations 
capture relatively well the effect of dissipation on the 
mutual diffusion coefBcient since the three approaches 
show a monotonic increase of D with decreasing a in all 
cases. On the other hand, at a more quantitative level, 
both first Sonine solutions overestimate slightly the pre- 
dictions of the second Sonine approach. In any case, the 
convergence of the Sonine expansion for this transport 
coefficient seems to be quite good, at least for not quite 
extreme values of mass and/or diameter ratios. 

We consider now the pressure diffusion coefficient 
Dp{a). This is plotted in Fig. 5 for the same cases as 
in Fig. 4. In contrast to the case of the mutual diffusion 
coefficient, when the defect species is lighter than the ex- 
cess component, the dependence of Dp on the coefficient 
of restitution predicted by the first Sonine approximation 
{Dp increases with decreasing a) differs from the one ob- 
tained from the more refined second Sonine solution [Dp 
decreases with decreasing a). At a more quantitative 
level, the first Sonine approximations overestimate again 
the second Sonine results for both values of the mass ra- 
tio, being the differences between both Sonine solutions 
more pronounced when /i < 1. Thus, at a = 0.5, the dis- 
crepancies between the first and second Sonine approxi- 
mations are about 4 % for = 4 while they are about 63 
% for fi = 0.5. The dependence of the thermal diffusion 
coefficient * on a is shown in Fig. 6. Note that, in the 
elastic limit, the first Sonine approximation to D * van- 
ishes while the second Sonine approximation is in general 
different from zero. We observe that both Sonine results 
tend to approach cac;li other as the dissipation increases. 
In particular, the dependence of D * on the coefficient 
of restitution predicted by the the first and second So- 
nine approximations is very weak when /i < 1 (in fact it 
is practically zero) while the coefficient increases clearly 
with dissipation in the opposite case {fj, > 1). In com- 
parison with the results obtained for Dp, the convergence 
of the Sonine solution for D' is better than that of the 
pressure diffusion coefficient, specially for strong dissipa- 
tion. It must be noticed that the differences between the 
standard and modified first Sonine approximations'*-^ are 
quite small in the region of coUisional dissipation consid- 
ered. Although not shown here, similar conclusions can 
be drawn when one considers other values for the mass 
and size ratios. 

As said in the Introduction, the results derived in this 
paper extend previous studies^^'^^ on both Sonine ap- 
proximations on the diffusion coefficients in the tracer 
limit {xi — >■ 0). Thus, one of the goals here is to as- 
sess the effect of finite concentration on the ratios of the 
second and first Sonine approximations to the diffusion 
transport coefficients. Figures 7, 8, and 9 shows the ra- 
tios D[2]/D[l], Dp[2]/Dp[l] and D'[2]/D'[l], respectively, 
versus the concentration xi for u = 1, a = 0.8 and two 
(disparate) values of the mass ratio /z. It is quite appar- 



ent that the impact of composition on the above ratios 
is in general significant. While the ratio I?[2]/£'[l] has 
a non monotonic dependence of xi, the corresponding 
ratios for the pressure and thermal diffusion coefficients 
exhibit a monotonic dependence with xi. Moreover, the 
second Sonine approximation to the diffusion coefficients 
differs clearly from its first Sonine approximation, spe- 
cially in the case of the thermal diffusion coefficient. 



VI. AN APPLICATION: THERMAL DIFFUSION 
SEGREGATION 

The knowledge of the diffusion transport coefficients 
allows us to examine thermal diffusion processes in binary 
mixtures. This is one of the most interesting problems 
appearing in multicomponcnt mixtures and it has been 
widely analyzed by some of the authors of the present 
paper^^'"^^'^^"'*^ in the tracer limit by considering the first 
and second Sonine approximations. The goal here is to 
assess the effect of the second Sonine solution on thermal 
diffusion segregation in the case of arbitrary concentra- 
tion Xi. 

We consider a granular binary mixture enclosed be- 
tween two plates at different temperatures. In a non- 
convecting steady state (u = 0) with gradients only 
along the orthogonal direction to the plates (z axis), 
the amount of segregation parallel to the thermal gradi- 
ent may be characterized by the thermal diffusion factor 
A. This quantity measures the separation of components 
caused by the temperature gradient. The factor A is de- 
fined as^'''^^ 



, dlnT d , /ni 
-A — - — = — in — 
dz dz V n2 



(51) 



Let us assume henceforth that cti > C72 and that the bot- 
tom plate is hotter than the top plate {dzT < 0). In this 
case and assuming that A is constant over the relevant 
ranges of temperature and composition, when A > the 
larger particles 1 tend to rise with respect to the smaller 
particles 2 (i.e., dz{ni/n2) > 0). In the opposite case, 
when A < the larger particles fall with respect to the 
smaller particles (i.e., dz{ni/n2) < 0). Although gravity 
is absent in our description, the former situation (A > 0) 
will be referred here to as the Brazil-nut effect (BNE) 
while the latter (A < 0) will be called as the reverse 
Brazil-nut effect (RBNE). 

Since no shearing fiows are present, the pressure ten- 
sor Pij = pSij and so the momentum balance equation 
(8) yields simply dzP = 0. Moreover, according to Eq. 
(7), ji.z = in the steady state. In the Navier-Stokes 
hydrodynamic order, ji^z is given by Eq. (1) so that the 
condition ji^^ = (along with dzP = 0) leads to the 
relation 



dlnxi 
dz 



P 



'2 D' d\nT 



mim2ni D dz 



(52) 
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The form of A can be easily obtained from Eqs. (51) and 
(52) and the result is 



A 



np 



D' 



1 D 



mim2'nin2 D X1X2 D* 



(53) 



where use has been made of the reduced expressions (42) 
for the mutual and thermal diffusion coefficients, respec- 
tively. Since the mutual diffusion coefficient D must be 
positive, the sign of A is determined by the sign of the re- 
duced coefficient D * . Consequently, the condition A = 
(which provides the criterion for the BNE/RBNE transi- 
tion) implies simply 

D* = 0. (54) 
According to Eqs. (45) and (46), the first Sonine approx- 
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FIG. 10: (color online) Plot of the thermal diffusion factor 
A[2] obtained from the second Sonine approximation as a 
function of the diameter ratio gi/g2 for an ordinary binary 
mixture (aij = 1) of hard spheres when both species have 
the same mass density (mi/m2 — (ai/o"2)^)- Three different 
values of the mole fraction are considered: (a) xi = 0.2, (b) 
xi = 0.5, and (c) x\ — 0.8. 

imation to Eq. (54) yields the criterion 

(l-l]=0. 

(55) 

In the elastic limit [aij = 1), C* = and so, A[l] = 
in the first Sonine approximation. However, away from 
the dilute gas limit, A[l] is not zero'^^''*'^ and segregation 
appears for ordinary mixtures. In the case of granular 
mixtures (a^j 7^ 1), the solution to Eq. (55) is simply'^^ 



mi 
m2 



T2 



(56) 



Note that if one assumes energy equipartition (Ti = T2), 
then segregation is only predicted for particles that dif- 
fer in mass, no matter what their diameters may be. It 
must be emphasized that the criterion (56) compares well 
with molecular dynamics simulations'*'* carried out in the 
tracer limit {xi — 0). 



The second Sonine approximation to Eq. (54) leads 
to a much more intricate criterion than Eq. (56). In 
particular, the results show that A [2] ^0 even for elastic 
collisions (a^j = 1).'^* This is illustrated in Fig. 10 for 
a binary mixture of hard spheres {d — 3) constituted by 
particles of the same mass density (mi/m2 — ((Ti/iT2)^)- 
In this case, A[2] is always positive and so, the larger 
particles tend to move towards the cold plate (BNE). 
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FIG. 11: (color online) Plot of the thermal diffusion factor 
A[2] as a function of the mole fraction xi for mi = m2, o"i = 
CT2 and different values of the coefficients of restitution: (a) 
ciii = a22 = 0.5, ai2 = 0.9, and (b) an = 0.8, 022 = 0.9, 
Qfi2 = 0.7. The solid lines correspond to the second Sonine 
approximation while the dashed lines refer to the first Sonine 
approximation. 

Another interesting limit case corresponds to the sit- 
uation in which segregation is only induced by inelas- 
ticity, namely, when one considers a binary mixture 
whose constituents differ only by their respective coef- 
ficients of restitution. This situation has been theoreti- 
cally studied** from the Boltzmann equation and it has 
been also confirmed'*^ by molecular dynamics simulations 
of two-dimensional binary mixtures. In order to analyze 
this effect. Fig. 11 shows the thermal diffusion factor A 
versus the mole fraction xi when mi = m2, <ti — 02 
and different values of the coefficients of restitution. As 
expected,*® we observe that segregation can occur due to 
inelasticity alone. It is quite apparent that in both sys- 
tems there is a change in the sign of A at a given critical 
value x\^c of the composition x\. Although the form of 
A differs in the first and second Sonine approximations, 
the value X\^c for each mixture is (practically) the same 
in both Sonine predictions. In the case (a), x\^c = 0.5 
due to symmetry considerations. 

Apart from the above limit situations, the dependence 
of A [2] on the parameter space is quite intricate. To as- 
sess the impact of the Sonine approximation on the form 
of the phase diagrams delineating the regimes between 
BNE and RBNE in the {<Ji/a2, mi/m2)-plane, let us con- 
sider again the physical case of hard spheres (d = 3) with 
a common coefficient of restitution {aij = a). Figure 12 
shows phase diagrams obtained from the first and sec- 
ond Sonine approximations for a = 0.8 and two values 
of the composition xi. Although the first Sonine approx- 
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FIG. 12: BNE/RBNE phase diagram for inelastic hard spheres at a — 0.8 and two different values of the mole fraction 
xi: a::=0.1 (panel (a)) and a;=0.5 (panel (b)). Points above the curves correspond to A > (BNE) while points below the 
curves correspond to A < (RENE). The dashed and solid lines are the results obtained from the first and second Sonine 
approximations, respectively. 
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FIG. 13: (color online) BNE/RBNE phase diagram for inelas- 
tic hard spheres with x\ = 0.7 and three different values of 
the (common) coefficient of restitution a. Points above the 
curves correspond to A > (BNE) while points below the 
curves correspond to A < (RBNE). 



over, in contrast to what happens in the dense case, at 
a given value of the concentration, the transition from 
BNE to RBNE may occur following two paths: i) along 
constant mass ratio mi/m2 with decreasing diameter ra- 
tio (Ti/tT2, and ii) along constant diameter ratio a\lu2 
with decreasing mass ratio m,\lm,i. 

Next, we study the impact of inelasticity on thermal 
diffusion segregation by considering only the most accu- 
rate theory (the second Sonine approximation). Figure 
13 shows the phase diagram for x\ — 0.7 and three val- 
ues of a (a =0.9, 0.7 and 0.5). The results show that 
the main effect of collisional dissipation is to reduce the 
size of the BNE region. This contrasts again with the 
results obtained from the first Sonine approximation for 
the Enskog equation (see, for instance. Fig. 5 of Ref. 43). 
The influence of dissipation on the BNE/RBNE phase 
diagram is much more significant for quite strong values 
of a (say for instance, a — 0.5) since the lines delineating 
the regimes between BNE and RBNE for a — 0.9 and 0.7 
are quite similar (at least in the region of values of the 
diameter ratio explored). 



imation reproduces qualitatively well the trends of the 
phase diagram, it clearly overestimates the predictions of 
the second Sonine approximation, specially at large size 
ratios and small mass ratios. The differences between 
both Sonine approaches becomes more important as the 
concentration of the larger particles decreases. These re- 
sults suggest that the Sonine expansion for the thermal 
diffusion factor exhibits a very poor convergence for di- 
lute gases in the tracer limit (x\ — 0). This conclusion 
agrees with previous results^^ obtained for moderately 
dense gases. Regarding the influence of the concentra- 
tion of the mixture x\ on phase diagrams we observe that 
the BNE region is reduced as x\ increases. On the other 
hand, this effect is less significant than for dense binary 
mixtures (see, for instance. Fig. 7 of Ref. 43). More- 



VII. SUMMARY AND DISCUSSION 

In this paper we have determined the mass flux j^^'' of a 
granular binary mixture at low-density. The results have 
been obtained by solving the inelastic Boltzmann equa- 
tion by means of the Chapman-Enskog method at the 
Navier-Stokes order. Three diffusion coefficients char- 
acterize the mass flux in the Navier-Stokes regime: the 
mutual diffusion coefficient D (that couples j^^"'^' with the 
concentration gradient Vxi), the pressure diffusion co- 
efficient Z?p (that couples j^^'' with the pressure gradi- 
ent Vp), and the thermal diffusion coefficient D' (that 

couples j^^'' with the temperature gradient VT). On the 
other hand, as for elastic collisions,^ the above coefficients 
[see Eqs. (16)-(18)] are defined in terms of quantities A.\^ 
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Bi, and Ci which are the solutions of a set of hnear inte- 
gral equations [see Eqs. (19)-(24)]. Given that the above 
quantities cannot be exactly obtained, they are approxi- 
mated by a truncated Sonine polynomial expansion. This 
allows us to obtain explicit forms for the diffusion trans- 
port coefficients in terms of the coefficients of restitution 
and the parameters of the mixture (relative masses, di- 
ameters and concentration). In this paper, we have de- 
termined D, Dp and D' by considering two polynomials 
in the Sonine polynomial expansion [see Eqs. (31)-(36)]. 
This approximation is usually referred to as the second 
Sonine approximation. Our present study complements 
and extends previous works on diffusion transport coef- 
ficients carried out in the tracer limit. ^^'^^ 

As mentioned in the Introduction, previous results^-'^ 
derived many years ago for ordinary mixtures (aij — 1) 
have clearly shown the reliability of the second Sonine 
approximation for the mutual and thermal diffusion co- 
efficients for a wide range of values of masses and sizes. 
These results have mainly encouraged the present work 
since the studies of the impact of the Sonine approxima- 
tion on the Navier-Stokes transport coefficients are very 
scarce in the case of granular mixtures. On the other 
hand, given the technical difficulties involved in the eval- 
uation of the second Sonine corrections to the transport 
coefficients, we have focussed on our efforts in the case 
of diffusion coefficients which are related to the lowest 
velocity moment (the mass flux) of the flrst-order distri- 
bution functions f^^"^- 

In order to gauge the accuracy of the second Sonine 
approximation, we have compared our theoretical pre- 
dictions for the mutual diffusion coefficient D with nu- 
merical solutions of the Boltzmann equation by means 
of the DSMC.^° Two situations have been considered: 
the self-diffusion problem (namely, when both species are 
mechanically equivalent) and the tracer limit (namely, 
when the concentration of one of the species is negligible). 
These are perhaps the two most simple situations where 
the coefficient D can be measured from the mean square 
displacement of a tracer particle immersed in a granular 
gas under HCS. The simulations performed here consider 
more general situations than those analyzed in previous 
works^^'^^'^^ where it was assumed that a22 = ol\2- As 
in previous studies, the present comparison shows again 
that in general the second Sonine approximation to D 
improves significantly the prediction of the first Sonine 
approximation, especially for high degree of dissipation 
and/or extreme mass or size ratios. 

The knowledge of the diffusion transport coefficients 
allows us to study segregation by thermal diffusion. This 
is a quite interesting problem not only from a fundamen- 
tal point of view but also from a more practical point 
of view. In a steady state without shearing flows, the 
thermal diffusion factor [defined by Eq. (51)] provides a 
measure of the separation or segregation induced by the 
presence of a temperature gradient. As Fig. 12 
clearly shows, the influence of the Sonine approximation 
over the form of the phase diagrams is quite important. 



especially for large values of the size ratio. In particu- 
lar, at a given value of cri/o"2, the flrst Sonine approach 
clearly overestimates the RBNE region. We expect the 
segregation criteria obtained here by using the second 
Sonine approximation can be tested against DSMC re- 
sults, molecular dynamics (MD) simulations, and even- 
tually experiments in real problems. We are currently 
working on DSMC and MD simulations adapted to the 
problem of segregation. 

One of the main limitations of the present study is 
its restriction to dilute gases. The Navier-Stokes trans- 
port coefficients for granular mixtures at moderate densi- 
ties have been recently obtained from the Enskog kinetic 
equation.^* These calculations have been done by consid- 
ering the first Sonine approximation. The extension of 
the results derived in this paper to higher densities is one 
of our possible future works. This extension will allow us 
to compare our theoretical results (based on the second 
Sonine approximation) with MD simulations performed 
at finite densities. 
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Appendix A: First and second Sonine 
approximations to the mass flux 

In this Appendix we determine the first and second So- 
nine approximations to the diffusion coefficients _D, Dp, 
and D' . These coefficients are defined by Eqs. (16)-(18) 
where the functions Ai, Bi, and C, are given by Eqs. 
(31) (36). Our goal is to evaluate the nine independent 
Sonine coefficients 

{ai,i; ''I,!; ci,i; ai,2; 012,2; &i,2; ''2,2; ci,2; 02,2}- (Al) 

The three first coefficients (oi^i, 61^1, and ci,i) are di- 
rectly related to the diffusion coefficients D, Dp, and D', 
respectively. 

Substitution of Eqs. (31)-(36) into the integral equa- 
tions (19), (21), and (23) gives 

_^(0) (T^T + pdp) [ai,iVi + ai,2Si(Vi)] /i,m 

-t-ai,i [£i/i,mVi - S-fMif2,MV2] 

+ai,2'Ci/i,MSi(Vi) + a2,2A^l/2,MS2(V2) 

= Ai2, (A2) 
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_^(0) (TdT+pdp) - 2C(°)] [61,1 Vi + 61,281 (Vi)] /i,M 

+61,1 [£i/i,mVi - SjMlf2,M^-2] 
+6i,2A/i,mSi(Vi) + 62,2A^i/2,mS2(V2) 

= B12, (A3) 



Bi2=Bi + ^-/i,m(ci,iV + ci,2Si), (A6) 



[ci,lVi+Ci,2Si(Vi)]/i,M 



+Cl,i [£i/i,mVi - S-fMlf2,M^2] 
+C1,2-Ci/i,mSi(Vi) + C2,2A1i/2,mS2(V2) 
= C12, 



(A4) 



where 
A12 



V 9X1 Jp,T 

+ T(ci,iV + ci,2Si)], 



/i,Mb(&i,iV + 61,281 



Ci2 = Ci-^/i,M (61,1V + 61,281). (A7) 



Here, Ai, Bi, and Ci arc given by Eqs. (27) (29), re- 
spectively. The corresponding counterparts of Eqs. (A2)- 
(A4) can be obtained from them by just making the 
change 1 o 2. Next, we multiply Eqs. (A2)-(A3) by 
(A5) TOiVi and integrates over the velocity. The result is 



-C(°) {TOt + pdp) + V niTiai^i + mTi (riiai,2 + ^202,2) = - ( ^"1^1) + (^) "1^1 ^P^^'^ + ^^^i) ' 

(A8) 



^(0) (T^T + pdp) - 2^' + V niTih^i + niT^ (rii6i,2 + ^262,2) = -'-^ ( 1 



>(o) 



niTi / minT\ TC,^^'^ 



-niTici,i, (A9) 



pCio) 

niTiCi,i + mTi (tii6i,2 + ri262,2) = — ^niTi6i,i. 



(AlO) 



Here, we have introduced the collision frequencies 

^ dVimiVi • [/:i/i,mVi - SjMif2,M^2] 



result is 



dniTi 
dn 



^ JdVimiVi ■ (ji2[vi|/i,MVi,/i°^] 

^7-/i2[vi|/1°\/2,mV2]), (All) 



InniTi 



+ TiiOi,2 + Ti2a2,2 



{ph,i + Tci,i) , 
(A14) 



(^i' - + '^11^1, 2 + ^1262,2 

1 / minT\ rC(°' 
P \ pTi J P 



(A15) 



Tii = J dv.miVi ■ A (/i.MSi) , (A12) 

Tij = j dvimjVi • Mi {fj,M^j) , i^j- (A13) 

^From dimensional analysis, niTiai,i ~ T^/^, niTi6i,i <^ 
T^f^ /p, and niTiCi,i ^ T~^/^. Thus, the temperature If only the first Sonine approximation is considered 
derivatives can be performed in Eqs. (A8)-(A10) and the (which means ai,2 = 6j,2 = Ci,2 = 0), the solution to Eqs. 



Z^Ci,i + riiCi,2 + Ti2C2,2 



pC 



2T 



-61,1. (A16) 
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(A14)-(A16) is 

Ol,l[l] = - 



d 

dxi J J, 



ci,i[l] = -^^'i,i[l], 



(A17) 

-1 

; 

(A18) 
(A19) 



Here, ai,2[l]. '^»,2[1], and Ci,2[l] denotes the first Sonine 
approximation to ai^2, bi,2, and Ci,2, respectively. From 
Eqs. (A17)-(A19) one gets the first Sonine expressions 
(44)-(46) for D, Dp and D' , respectively. 

To close the problem, we multiply now Eqs. (A2)-(A4) 
by Si(Vi) and integrates over the velocity. Following 
identical mathematical steps as before and after some 
algebra one gets 

^Z^ll - ^C^°^^ «1,2 + !^12a2,2 - (~^^^ 



+Tci,2) 



Here, A12 ~ An — 6^X12, where 6 = x\/x2, X2 = 1 — Xi, 
7 = T1/T2 and 



d{d + 2) mTf 



dw^Si ■ A (/i,MVi) , (A23) 



A, 



(A24) 

In addition, in Eqs. (A20)-(A22), we have introduced the 
collision frequencies 




d{d+2) n^Tf 



dviSi ■ Mi (/j.mSj) , 



(A25) 
(A26) 

The corresponding integral equations verifying the re- 
maining coefficients a2.2, &2.2, and C2.2 can be obtained 
from Eqs. (A20) (A22), respectively, by interchanging 
1 O 2. Note that upon writing Eqs. (A20)-(A22) we 

have neglected the non-Gaussian corrections to /^''^ . 

Equations (A14)-(A16) along with Eqs. (A20)-(A22) 

can be written in a more compact form by using ma- 
trix notation. For the sake of convenience, let us intro- 
duce the dimensionless coefficients a* = voai i, fe* ^ = 
piyobiA, c* 1 = Ti^oci,!, a*2 = Ti'oai,2, 6*2 = pTi'oha^ 
and c*r, = T^j/oCi,2, where is defined by Eq. (43). Let 



(A20) us introduce the column matrix X by 



'2,2 Ci 2 

P 



V(0) 

- ( ^ - A12 ) 61,1 = 0, 



{'^iji! ^1,1) ^1,1; c*i,2> '^2,2! ^1,2! ^2,2> ^1,2! ^2,2}- (A27) 

Therefore, according to Eqs. (A14)-(A16) and (A20)- 
(A21) (A22) , the coupled set of nine equations for the unknowns 
can be rewritten in matrix form as 



(^11 - c 



(0) 



Cl,2 + J^12C2,2 



2T 



Cio) \ 1 

— A12 ) Ci 1 = — — — - . (A22) where the square matrix ft is 

Ti / ' TTi 



(A28) 



f2 = f2(°)+f|(i). 



(A29) 
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(V* 



2 ^ 








T* 


/no 
±z 











\ 





2 ^ 











11 


T* 

' 12 














V* 














Tn 













11 2 ^ 


12 























Zi 

































1^12 

























^^22 2'' 





























^li - c 


^12 























1^21 





( ^ 


_(0C\ 
















-c 











C/2 











He 

Wl2 














H= 

^21 

















'^12 














'^21 

















^t2 








V 





^21 









The column matrix Y is given by 



ac] 





(OC] 























-c 














-c* 


c/2 














c/2 





/ 



Xi7i 
-d- 



\ dxi 



1 

2f 'f \ Uxi 

1 (djl 
2^ \dx 






ji{i+d)_\ 

■yi{l+IJ.6) J 




-7i_ 
-72 



(A31) 



(A32) 



In the above equations, we have introduced the reduced 
quantities v* = v/vq, t^j = Tij/TvQ, = Uij/uQ, and 



C* ( C 

ujl2 = ^*i2- — , ^^21 = -^7 - — 

7i V 72 

The solution to Eq. (A28) is 



A* 



(A33) 



(A34) 



^From this relation one gets the second Sonine correc- 
tions to the coefficients an, 6ii, and cn. 



Appendix B: Reduced collision frequencies and 
cooling rates 

In this Appendix we provide the exphcit expressions 
of the (reduced) collision frequencies needed to evaluate 



-D[2], -Dp [2], and D'[2]. As said in the main text, to eval- 
uate them we take the local Maxwellian approximations 
(30) for the zeroth-order distributions ff'\ These col- 
lision frequencies have been already evaluated in the d 
dimensional case.^^'^^ They are given by 



27rl-i 

^ = -^7dV^^ + "i2) 



rfr(f) 



h+02 
$102 



1/2 



(a;2M2i +X11112) , 
(Bl) 



2:2(1 + ai2)- 



61 



3/2 



7r2 



T12 



1:2(1 + 012)7' 



-I <'(^l+g2)-^/^ 



(B2) 



, (B3) 
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'12 



r(f) d 



0'12 




(^1 + 02)"'/' (0102)"^/' (a;2A - -yxiB) , 



(B4) 



7r2 



'12 



1 ^,^1(1 + ^,12) 



r(f) d(rf + 2) Vi2 



Oi 



3/2 



F + (d + 2) 



01 + 02 



B 



(B6) 



TT^-^ 8 / o^y 

r(f) d(rf + 2) Vai2y 

^ + Y^(^^ + 8)(l-aii) 

7r('«-i)/2 1 
H — „ . . „>, a:2M2i(l + 0:12) 



a;i(2ei)-i/2(l + an) 



r(f) rf(rf + 2)' 



01 



\ 3/2 



02(0 



1 + 02); 



In the above equations, inj = mi/{mi + mj), and 



,=1 



(B7) 



In addition, the quantities A, B, E, and F are given, 
(B5) respectively, as 



J 



A = {d + 2)(2/3i2 + 02) + /i2i(0i + 02) {{d + 2)(1 - ai2) - [(11 + d)ai2 -M- 7]/3i20r'} 

+3(^ + 3)^2^01-1 + 2/xii (^2af2 - ^"12 + <i + l) 0r'(0i + 02)' - (rf + 2)020r'(0i + 02), 

B = {d + 2)(2/3i2 - 0i) + Ai2i(0i + 02) {{d + 2)(1 - 012) + [(11 + d)ai2 - 5rf - 7]/3i202~i} 
-3(d + i)Pl2e2^ - (2al2 - ^-^ai2 + d+l] e^^Oi + 62)^ + {d + 2)(0i + 62), 



E = 2f4^9^^ei + 62^ (^2al^ - ^-^ai2 +d+l^[{d + 2)9i + {d + 5)62] 

-M2i(0i + 02) {/3i20r'[(d + 2)01 + {d + 5)02][{n + d)ai2 -5d-7] 

-020rM2O + d{15 - 7ai2) + d^{l - aia) - 28^12] - (rf + 2)^(1 - 012)} 

+3{d + 3)l3f29i^[id + 2)^1 + {d + 5)02] + 2/3i20rM(c^ + 2)'0i + (24 + lid + d'^)92] 

+{d + 2)020r' [(«^ + 8)01 + + 3)02] -{d + 2)(0i + 02)0r'02[(d + 2)01 + {d + 3)02], 



F = 2iil,9^^9, + 02)2 (^2al2 - ^c.12 + + 1^ [(d + 5)0i + (d + 2)^2] 

-A'2i(0i + 02) {Pi202 ^[{d + 5)01 + {d + 2)02][(11 + rf)ai2 - 5d - 7] 

+0i02"i[2O + rf(15 - 7ai2) + rf2(l - 0,12) - 28ai2] + (rf + 2)^(1 - 0,13)} 

+3(0! + 3)/3f202"' [{d + 5)01 + (d + 2)612] - 2f3i202 ^ [(24 + llrf + d^)0i + {d + 2)^02] 

+ id + 2)9i0^^[{d + 3)01 + {d + 8)02] -{d + 2){9i + 02)92 ^[{d + 3)0i + {d+ 2)02]. 



(B8) 



(B9) 



(BIO) 



(Bll) 



Here, /3i2 = /Lti202 — At2i0i- From Eqs. (B2)-(B6), one by interchanging 1 o 2. 
easily gets the expressions for T22, T21, Aji, 1^22 and 1/21 
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Finally, the temperature ratio 7 is determined from the 
condition^^ 



C* 
1 — '=2 — S ) 

where the dimensionless cooling rate is 

\/27r 



(B12) 



4:TT~ 



1/2 



1- -^(1+Q!12) 



IC2 

^1+02 



(1 + ai2) 



(B13) 



The expression of Q can be obtained form the change 
102. Once the temperature ratio 7 is known, the 
partial temperature ratios ji = Ti/T {i = 1,2) can be 
expressed in terms of the (global) temperature as 



71 



l + a;i(7- 1)' 



72 = 



l + a;i(7-l)' 



(B14) 
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